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SUMMARY 


An  experimental  and  analytical  study  on  the  effect  of  thin  ceramic  coatings  on  soot 
and  NOx  emissions  and  performance  of  a  direct  injection  dksei  engine  was  conducted 
jointly  between  North  Carolina  A&T  State  University  and  North  Carolina  State 
University. 

Performance  and  emissions  data  were  gathered  on  a  normally  aspirated  Ricardo 
Hydra  single  cylinder  DI  engine  with  various  combinations  of  ceramic  coatings  instaUed. 
Thin  ceramic  thermal  barrier  coatings  were  applied  to  the  piston  crown  and  bowl,  the 
head  and  valves,  and  the  cylinder  liner.  The  coated  piston  and  head  were  run  singly 
and  in  combination  with  the  cylinder  liner  to  investigate  the  effects  of  these  different 
coated  surfaces  on  emissions  and  performance  for  two  different  pure  hydrocarbon  fuels, 
hexadecane  and  dodecane.  Coating  the  piston  crown  alone  results  in  generally  lower 
cylinder  pressure,  lower  brake  specific  fuel  consumption  and  lower  NOx  emission 
compared  to  the  baseline  engine.  Soot  emission  is  typically  increased  below  2000  RPM 
and  decreased  above  2000  RPM.  Coating  the  head  alone  reduces  cylinder  pressure,  but 
generally  increases  specific  fuel  consumption  and  NOx  and  soot  emission. 

The  analytical  portion  of  the  study  involved  modifications  to  the  KIVA-H  code  and 
its  use  to  model  the  Hydra  engine  with  the  thermal  coatings.  Modifications  to  the  code 
include  incorporation  of  an  eddy-break*up  combustion  model,  to  replace  the  standard 
Arrhenius  single  reaction  model  A  time  dependent  combustion  surface  temperature  was 
also  incorporated  to  simulate  the  effects  of  thermal  barrier  coatings  on  blinder 
temperature,  pressure  and  NOx  production.  A  soot  model  was  added  to  the  code, 
following  the  work  of  Magnussen  and  Hjertager.  The  EBU  model  gives  better  results  for 
the  diffusional  portion  of  the  combustion  process,  but  faib  to  adequately  model  the 
premixed  combustion,  typically  resulting  in  a  lower  predicted  peak  cylinder  pressure  (and 
temperature)  th—  predicted  by  the  single  reaction  model  and  shown  by  experiments. 
Consequently,  the  EBU  model  also  nnder>predicts  NO  emission  to  a  greater  degree  than 
the  SR  model  The  KEVA-H  modeling  has  led  to  an  understanding  of  the  effect  of  coating 
the  piston  on  NO  production.  The  hotter  piston  crown  warms  the  intake  air,  shortening 
ignition  delay  and  decreasing  the  ratio  of  premixed  to  diffhsion  combustion,  ultimately 
resulting  in  lower  peak  cylinder  temperature  and  reduced  NO.  The  lOVA-H  results 
agree  reasonably  well  with  the  experimental  data  for  cylinder  pressure  and  NO  and  soot 
emission. 


i 


LIST  OF  PUBUCATIONS 


Technical  Ptmcrs  PnbUshcd 

Afiiy,  E.Mm  Klctt  D.E.,  SiB^  H.  and  Wei,  X.«  "The  Effect  of  Inwhitfaig  the  Piston  Crown 
on  the  Performance  and  NO  Emissions  of  a  DI  Diesel  Engine'',  Proceedings  of  the  Ei^th 
International  Conference  for  Mechanical  Power  Engineering,  Alexandria,  Egypt,  April, 
1993. 

Saad,  M.,  Aliiy,  E.M.,  Klett  D.E.  and  Singh,  H.,  "Modding  of  Combustion  and  NO 
Emission  in  a  Direct  Injcctkm  Dksd  Engine",  Proceedings  of  the  Ei^th  International 
Conference  for  Mechanical  Power  Enghtccring,  Alexandria,  Egsrpt,  April,  1993. 

Saad,  M.,  Afify,  E.M.  and  Kktt,  D.E.,  "MnhidimcBsioBal  Modding  of  Combustion  and 
NO  Emission  of  the  Ricardo  Hydra  Direct  InJectioB  Diesel  Engine:  KIVA-n  Versus 
Experiment",  Proceedings  of  the  ASME  1C  Engine  DivisioB  Fall  Technical  Conference, 
Morgantown,  WV,  September,  1993. 


Papers  in  PreparatioB 

"Effect  of  Sdective  InsnIatioB  on  the  CombnstioB,  Performance  and  NO  Emission  of  DI 
Diesel  Engines" 

"Modeling  of  Soot  and  NO  Emimjeni  in  a  DI  Diesd  Engine:  KIVA-II  Versus 
Experiments" 


Theses 

Saad,  M.T.,  "Soot  and  NOx  Emisdeni  and  Combustion  Characteristics  of  Direct  Injection 
Diesel  Engine",  PhD  ThesiB,  Mechanical  Engineering  Department,  North  Carolina  State 
University,  December,  1993. 

Vakil,  M.,  "Experimentd  Investigation  of  Effect  of  Ceramic  Coatings  on  Performance 
and  Emissions  of  a  DI  Diesel  Engine",  MS  Thesis,  Mechanical  EnginccriBg  Department, 
North  Carolina  A&T  State  University,  (In  Progress) 


ii 


i 


LIST  OF  PERSONNEL  AND  DEGREES  EARNED 


David  E.  Klctt,  Project  DIrctor,  Profenor  of  Mcchaakal  Eagtaccriag,  North  Caroliaa  A&T  State  University 
E.  M.  AfUy,  Co-Investigator,  Professor  of  Mechanicai  Engineering,  North  Carolina  State  University 
Messiha  T.  Saad,  Gradnate  Research  Assistant,  PhD,  North  Carolina  State  University,  12/93 
Xu  Wei,  Gradnate  Research  Assistant,  PhD,  North  Carolina  State  University,  Expected  S/95 
M.  Vakil,  Gradnate  Research  Assistant,  MS,  North  Carolina  AAT  State  University,  Expected  5/94 
Nigel  Brooks,  Undergraduate  Research  Assistant,  BS,  North  Carolina  AAT  State  University,  Expected  12/94 
Martin  Ashley,  Undergraduate  Research  Assistant,  BS,  North  Carolina  AAT  State  University,  5/92 
Raphael  Jones,  Undergraduate  Research  Assistant,  North  Carolina  AAT  State  Univcfsity,  5/93 
Postoria,  Aguirre,  Undergraduate  Research  Assistant,  North  Carolina  AAT  State  University,  5/92 


iii 


ACKNOWLEDGMENT 


TUs  work  was  supported  Iqr  the  U.  S.  ^nqr  Rcscardi  Oflke  under  grant  number  DAAL4I3- 
8Mr4)l<l3.  D.  E.  Kktt  was  dm  project  director  and  principal  inrestigator  id  North  CaroUna  AAT 
State  Unhrersity.  E.M.Afify  was  the  principal  inrestigatMr  at  North  Carolina  State  Unirersity,  and 
M.  Saad  was  a  research  associate  on  the  project.  The  authors  would  like  to  express  their 
appreciation  to  Darid  Man  of  die  Army  Research  Office  for  his  assistance,  adrke  and 
encouragement  during  the  dnradon  of  dds  project. 


1.  INTRODUCTION 


1.1  PROBLEM  STATEMENT 

The  objective  of  this  work  was  to  study  the  effects  of  thin  ceramk  coatings  applied 
to  various  internal  engine  surfaces  on  the  performance  and  emissions  of  a  direct  injection 
(DI)  diesel  engine.  The  study  included  an  experimental  investigation  coupled  with 
computer  modeling  using  the  KIVA  II  code  to  gain  insight  into  the  combustion  effects 
caused  by  higher  surface  temperatures  resulting  from  the  use  of  ceramic  coatings.  The 
project  was  a  joint  effort  between  North  Carolina  A&T  State  University,  where  the 
experiments  were  carried  out,  and  North  Carolina  State  University,  where  the  bulk  of  the 
computer  modeling  was  performed. 

The  Army,  among  a  number  of  other  agencies  and  corporations,  has  maintained 
an  interest  in  the  use  of  ceramics  in  diesel  engines  for  a  number  of  years,  and  many 
studies  have  been  conducted  on  various  aspects  of  this  topic  [1-8).  This  study  focused  on 
the  use  of  KIVA  II  to  model  the  combustion  process  in  a  particular  DI  diesel  engine  with 
thin  ceramic  coatings  applied  to  the  piston  crown,  the  cylinder  head  inner  surface  and  the 
cylinder  liner.  Experiments  were  performed  on  a  single-cylinder  Ricardo  Hydra  DI  diesel 
research  engine  to  provide  a  data  base  for  comparison  with  the  modeled  results.  Of 
particular  interest  were  the  effects  of  the  coatings  on  brake  specific  fuel  consumption, 
NOx  emission  and  soot  emission.  Two  pure  hydrocarbon  fuels,  Dodecane  and 

Hexadecane  were  used  for  the  experiments  and  in  the  combustion  modeling. 

These  fuels  were  chosen  for  purposes  of  reproducibility  of  experimental  results  and 
because  the  properties  are  available  to  facilitate  the  modeling. 

1.2  BACKGROUND 

The  use  of  ceramic  coatings  on  combustion  chamber  surfaces  influences  the 
performance  and  nhanst  emissions  of  DI  diesel  engines  [1-4].  This  can  be  attributed  to 
changes  in  the  combustion  process  brought  about  by  the  modification  of  the  thermal 
boundary  conditions  due  to  the  lower  thermal  conductivity  of  ceramics  relative  to  metals. 
The  added  thermal  resistance  of  thin  low  conductivity  coatings  increases  the  surface 
temperature,  and  higher  surface  temperatures  affect  the  cowustion  process  in  a  variety 
of  ways,  depending  on  which  combustion  chamber  surfaces  are  coated. 
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To  better  understand  the  effect  of  insulating  die  combustion  chamber  walls  on  the 
performance  and  exhaust  emissions  of  the  engine,  a  brief  review  of  the  combustion 
process  in  DI  diesel  engines  is  first  presented. 

The  combustion  process  in  Dl  diesel  engines  involves  two  modes  of  combustion, 
usually  referred  as  the  premixed  and  the  diffusion  combustion  modes  [9].  Premixed 
combustion  occurs  early  in  the  process  when  fuel  which  has  evaporated  and  mixed  with 
air  during  the  ignition  delay  period  auto  ignites.  This  mode  is  accompanied  by  a  high 
rate  of  heat  release  which  produces  a  rapid  rate  of  pressure  rise.  When  the  premixed 
fiiel/air  mixture  is  depleted,  diffusion  combustion,  characterized  by  a  lower  rate  of  heat 
release,  takes  over  and  controls  the  remainder  of  the  combustion  process. 

Modification  of  the  engine  design,  or  variation  of  its  operating  conditions,  will 
affect  the  combustion  process  in  such  a  way  as  to  vary  the  ratio  of  the  time  duration  of 
the  two  modes,  tp/t^,  where  tp  represents  the  duration  of  the  premixed  mode  and  t^ 
represents  the  duration  of  the  diffusion  mode.  This  ratio  plays  an  important  role  in 
determining  the  effect  of  engine  parameters  and  other  factors  on  the  performance  and 
emissions  of  the  engine.  Factors  which  increase  the  ignition  delay  period  will  increase 
the  premixed  combustion  duration  while  decreasing  the  diffusion  combustion  duration, 
causing  an  increase  in  the  ratio  tp/t^.  This  may  lead  to  higher  peak  cylinder  pressure 
and  temperature  which  may  improve  thermal  efficiency  and  reduce  CO  and  unbumed 
hydrocarbon  (UHC)  emissions  at  the  expense  of  increasing  NOx  emissions  of  the  engine. 
Large  increases  in  tp/t^  will  cause  a  high  rate  of  pressure  rise  and  lead  to  objectionable 
diesel  knock.  Factors  such  as  engine  speed,  advanced  injection  timing  and  the  use  of  low 
cetane  fiiels  contribute  to  longer  ignition  delay  and  the  increase  in  tp/t^.  Reducing  the 
ignition  delay  period  causes  the  premixed  combustion  duration  to  decrease  while 
increasing  the  diffusion  combustion  duration,  Le.  reduces  tp/t^,  prolonging  the  combustion 
process.  Large  decreases  in  tp/t^  may  cause  the  diffusion  combustion  mode  to  dominate 
which  may  lead  to  loss  of  power,  decrease  in  thermal  efficiency  and  possible  deterioration 
of  engine  exhaust  emissions.  Preheating  the  intake  air  and  insulating  combustion 
chamber  surfaces  are  among  the  factors  which  can  decrease  ignition  delay  and  tp/t^. 

Several  experimental  studies  on  low  heat  rejection  (LHR)  DI  naturally  aspirated 
diesel  engines  (13>4,5,6]  showed  that  engine  performance  and  exhaust  emissions  suffered 
from  lengthening  of  the  combustion  process  which  is  indicative  of  an  extreme  reduction 
in  the  ratio  tp/t^. 
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1.3  PERTINENT  LITERATURE  REVIEW 


Few  investigators  have  studied  the  effect  of  insulating  various  combustion  chamber 
surfaces  in  DI  diesel  engines  and  the  results  have  been  contradictory  in  some  cases. 
Miyairiy  et.  al  [3]  studied  the  effect  of  selective  insulation  of  the  cylinder  head,  piston 
crown,  and  cylinder  liner,  using  thick  monolithic  ceramic  inserts,  on  the  performance  and 
emission  characteristics  of  a  single-cylinder,  normally-aspirated  DI  diesel  engine.  They 
showed  that  fuel  economy  and  NO  emissions  of  the  engine  were  improved  by  insulating 
the  cylinder  head  and  liner,  but  were  made  worse  by  insulating  the  piston  crown.  Part 
of  the  degradation  in  BSFC  with  the  insulated  piston  was  attributed  to  the  increased 
reciprocating  mass  due  to  the  heavy  monolithic  ceramic  piston  crown. 

Assanis,  et.  al  [2]  conducted  a  series  of  tests  on  a  supercharged  DI  diesel  engine 
with  and  without  puton  surface  insulation  to  determine  the  effect  of  ceramic  coating  the 
piston  crown  on  engine  performance  and  emissions.  In  their  study,  they  emphashced  the 
significance  of  the  heat  release  profile,  and  indicated  that  insulating  the  piston  with  a  thin 
coating  of  PSZ  resulted  in  better  engine  efficiency  and  reduced  emissions  over  the 
baseline  engine. 

Dickey  [I]  studied  the  effect  of  applying  thin  ceramic  coatings  to  all  combustion 
chamber  surfaces  in  a  supercharged  single-cylinder  Caterpillar  lY-540  DI  diesel  engine. 
The  results  showed  decreased  thermal  efficiency,  but  also  decreased  specific  NOx  and 
UHC  for  the  ceramic  coated  engine  relative  to  the  baseline  engine,  especially  at  higher 
loads. 


Daby,  et.  al  [10]  conducted  tests  on  a  single-cylinder,  uncooled,  pressurized-intake, 
DI  diesel  engine  equipped  with  thick  monolithic  ceramic  inserts  on  the  head  and  cylinder 
liner  (above  top  ring  traveO*  The  piston  was  a  low-heat-pass  articulated  design  with  a 
ferrous  metal  ti^  and  an  aluminum  skirt.  These  tests  were  conducted  to  investigate  the 
behaviour  of  altcmative  fuels  (low  cetane,  higher  density)  in  conjunction  with  an  uncooled 
engine.  Basline  tests  were  conducted  with  Phillips  diesel  control  fuel  D-2  in  the  uncooled 
engine  and  in  a  baseline  cooled  metal  engine.  The  baseline  engine,  however,  had  a 
different  stroke  and  compression  ratio  than  the  ceramic  engine,  and  direct  comparisons 
of  performance  and  emissions  for  the  two  engines  would  not  be  meaningfuL 


2.  DESCRIPTION  OF  EXPERIMENTS 


2.1  EXPERIMENTAL  APARATUS 

The  experimental  data  used  for  comparison  with  the  KIVA  results  were  obtained 
from  a  Ricardo  Hydra  single^ylinder  D1  diesd  research  engine.  Table  1  provides  the 
basic  specifications  of  the  engine  and  Figure  2.1  provides  a  cross-sectional  view  of  the 
combustion  chamber  showing  the  toroidal  bowl  piston.  Measured  data  included  load, 
speed,  crank  angle,  needle  lift,  cylinder  and  fuel  line  pressure,  fiiel  mass  flow  rate,  air 
volume  flow  rate,  intake  and  exhaust  gas  temperatures,  coolant  and  lube  oil  temperatures, 
exhaust  smoke  opacity  and  NO  and  NOx  emissions. 

Engine  loading  was  via  a  digitally-controlled  DC  motoring  dynamometer 
interconnected  to  the  utility  grid  through  a  KTK  power  control  and  signal  conditioning 
system.  Speed,  load  and  injection  timing  could  be  remotely  controlled  from  an 
instrumentation  panel  which  provides  readouts  of  these  parameters  plus  coolant, 
lubricant,  intake  and  exhaust  manifold  temperatures.  Additional  digital  thennoconples 
were  used  to  monitor  room  air  temperature,  NOx  sampling  line  temperature  and  air  filter 
intake  temperture.  A  sUng  p^chrometer  was  used  to  determine  combustion  air  humidity. 
The  engine  coolant  temperature  and  oil  temperature  could  be  individual^  controlled  with 
a  pair  of  thermostats  that  r^ulate  the  flow  of  laboratory  water  through  heat  exchangers 
that  provide  cooling  for  the  engine  fluids. 

Exhaust  smoke  opacity  was  measured  with  a  USPHS  Diesel  Smokemeter 
permanently  mounted  in  the  exhaust  pipe  75  mm  downstream  from  the  exhaust  port. 
The  smoke  meter  was  calibrated  against  a  set  of  neutral  density  geUtin  filters.  Oxides 
of  nitrogen  were  measured  with  a  Thermo  Environmental  Instruments  Model  lOAR 
chemiluminescent  NO-NOx  analyzer  in  conjunction  with  a  Model  800  heated  sample 
conditioning  unit  and  a  heated  sampling  line.  The  NOx  system  was  calibrated 
periodical^  against  2000  ppm  NO  in  nitrogen  calibration  gas. 

Intake  air  flow  was  measured  with  a  Meriam  Model  50MC2-2F  Laminar  Flow 
Element  equipped  with  a  digital  manometer  to  read  the  pressure  differential,  and  an 
inclined  water  manometer  for  calibration  purposes.  Fuel  mass  flow  rate  was  monitored 
with  an  AVL  Model  730  gravimetric  fuel  balance. 
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The  cylinder  pressure,  fuel  pressure,  needle  lift  and  crank  angle  were  recorded 
with  a  Data  Precision  DATA  6000  digital  wave  form  analyzer  at  a  data  rate  of  25  KHz. 
The  data  was  stored  on  diskettes  for  subsequent  computer  analysis.  Table  2  gives  the  test 
matrix  of  load,  speed  and  injection  timings  used  for  the  experiments. 

Ignition  delay  and  injection  duration  were  determined  from  the  digital  wave  form 
data.  The  25  KHz  data  rate  translates  to  40  microseconds  per  data  point,  which  then 
represents  the  uncertainty  for  both  the  ignition  delay  and  the  injection  duration. 


Figure  2.1 

Cross  Section  of  Ricardo  Hydra  DI  Diesel  Combustion  Chamber 
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TABLE  1  -  SPECIfICATIOI«S 


Swept  Volume 


Maximum  Speed 


Maximum  Power 


Maximum  Cylinder  Preuure 


Compreuion  Ratio 


Connecting  Rod  Length 


Squish  Height 


Valve  Timing: 


Intake  Opens 


Intake  Closes 


Exhaust  Opens 


Exhaust  Closes 


InJocHaa  Syalem 


Injector  Pump 


Noxzle  Opening  Pressure 


Injector 


Lift  Pump 


Manu&etaror 


Rating 


Max.  Speed 


120  bar 


0.82542  mm 


42  ATDC 


S8BBDC 


10  ATDC 


Micro  Bosch  sixe  A  type  EA 


4  holes  *  0.21  mm  dia.  *  155 


250  bar  »  25  Mpa 


KBEL  88  PV  1187 


McClure 


Shunt  wound  dc  wHh  separate  excitation 


30  KW  continuous 


100  rps 


KTK  type  6P4Q3D  converter  for  motoring  and 
regenerative  loading 
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2.2  CERAMIC  COATINGS 


The  piston,  cylinder  hend,  valves  and  cylinder  liner  were  sent  to  Adiabatics,  Inc. 
in  Columbus,  IN  for  coating.  The  piston  crown  and  bowl  received  a  0.25  mm  slurry- 
sprayed  coating  of  partially  stabilized  zirconia  (PSZ)  consisting  of  85  percent  partially 
calcium  stabilized  cubic  zirconia,  10  percent  tungsten  cobalt  chrome  powder,  and  5 
percent  chrome  oxide.  The  head  and  valves  were  coated  with  a  0.5  mm  thermal  barrier 
coating  incorporating  5  percent  hollow  alumina  spheres  in  a  slurry  of  65  percent  silica, 
15  percent  PSZ,  7  percent  tungsten  chrome  powder,  and  8  percent  chrome  oxide. 

The  cylinder  liner  was  bored  0.2  mm  over  the  entire  length  plus  an  additional  0.75 
mm  in  the  region  above  the  top  ring  reversal  (TRR).  The  region  above  the  TRR  was 
plazma  sprayed  with  ytria  stabilized  zirconia,  and  then  the  entire  length  of  the  liner  was 
given  a  0.2  mm  wear  coat  of  slurry  sprayed  PSZ.  The  coatings  were  sealed  and  densified 
with  a  chrome  oxide  based  drain  cast  slurry. 

2.3  TEST  MATRIX 

A  full  set  of  data  runs  were  conducted  with  each  of  the  two  fuels,  hexadecane  and 
dodecane,  and  for  each  of  four  different  en^e  builds  utilizing  different  combinations  of 
coated  components  including:  (1)  baseline  (no  coatings),  (2)  coated  piston  alone,  (3)  coated 
head  alone,  and  (4)  coated  piston,  head  and  liner  together.  A  full  set  of  runs  consisted 
of  the  test  matrix  shown  in  Table  2  involving  four  engine  speeds,  four  loads  and  four 
injection  timings,  for  a  total  of  64  runs  per  engine  build,  per  fuel  (512  total  runs).  After 
initially  investigating  the  effect  of  engine  coolant  temperature  on  NOx  production  and 
observing  the  expected  increase  in  NOx  with  increasing  temperature,  as  shown  in  Figure 
2,  it  was  decided  to  use  a  constant  coolant  temperature  of  353  K  and  constant  oil 
temperature  of  323  K.  A  sample  data  set  for  hexadecane  for  one  engine  speed  is  shown 
in  Table  3. 


1 _ _ 

INJECTION  TIMING  (DEG  BTDC) 

16 

18 

20 

22 

BMEP  (BAR) 

2.24 

3J5 

3.91 

4.47 

Speed  (RPM) 

1000 

1500 

2000 

2500  1 
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TABLES  EXPERIMENTAL  DATA 

FUEL:  BEXADECANE 

FUEL  DEF^TY:  0.773  gp/en’ 

ISeORPM 

RUN  NUMBER  | 

TEST  DATA 

1 

2 

3 

4 

5 

D 

7 

8 

LOAD  N-m 

8 

8 

8 

8 

12 

12 

12 

12 

INJ.  TIMING 

DEG-BTDC 

16 

18 

20 

22 

16 

18 

20 

22 

INTAKE  AIR  TEMP.  K 

310 

311 

312 

312 

312 

312 

312 

312 

WATER  TEMP.  K 

353 

353 

353 

353 

353 

353 

353 

353 

EXHAUST  TEMP  K 

499 

505 

515 

519 

545 

551 

560 

565 

SMOKE  LEVEL  (Amp.) 

48 

47 

47 

46 

45 

45 

43 

42 

SOOT  AT  STP  (g/itf ) 

.234 

.234 

.300 

.335 

.367 

.499 

.565 

.635 

SOOT  AT  EXHAUST 
(g/m^) 

.141 

.138 

.176 

.372 

.202 

.270 

.304 

.372 

1  NO,  (PPM) 

1300 

1550 

1750 

1900 

1900 

2150 

2300 

25. 

NO  (PPM) 

1125 

1250 

1450 

1600 

1650 

1850 

2000 

2150 

FUEL  INJ.  g/min 

8.2 

8 

8.3 

8.6 

10 

10.2 

10.5 

10.9 

IGNIT.  DELAY  (msec) 

0.68 

0.72 

0.8 

0.84 

0.56 

0.64 

0.72 

0.76 

INJ.  DURATION  (msec) 

o.% 

1.04 

1.08 

1.08 

1.12 

1.2 

1.2 

INJ.  DURATION  (DEG) 

8.64 

9.0 

9.36 

9.72 

9.72 

10.1 

10.8 

10.8 

FUELLING  RATE  *10-* 

(g/inj) 

1092 

1066 

1106 

1146 

1332 

1359 

1399 

1452 

IGN.  DELAY  DEG 

6.12 

6.48 

7.2 

7.56 

5.04 

5.76 

6.48 

6.84 

BMEP  (bar) 

2.24 

2.24 

2.24 

2.24 

3.35 

3.35 

3.35 

3.35 

POWER  OUTPUT  W 

1257 

1257 

1257 

1257 

1885 

1885 

1885 

1885 

BSFCg/kwh 

391.5 

382.0 

396.3 

410.6 

318.3 

324.7 

334.2 

347.0 

AIR  FLOW  RATE 
(kg/b) 

21.25 

21.25 

21.25 

21.25 

21.25 

21.25 

21.25 

21.25 

AIR  FUEL  RATIO 

42.65 

44.28 

42.68 

41.19 

35.42 

34.73 

33.73 

32.50 

VOL.  EFFICIENCY  % 

91.38 

91.38 

91.38 

91.38 

91.38 

91.38 

91.38 

91.38 

ROOM  TEMP.  F 

72 

72 

72 

72 

72 

72 

72 

72I 
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TABLES 


#.773gM/Ga* 


:  I  ;J  D( ' 


10 

11 

12 

13 

14 

IS 

16 

LOAD  N-m 

14 

14 

14 

14 

16 

16 

16 

16 

INJ.  TIMING 

DEG-BTDC 

16 

18 

20 

22 

16 

18 

20 

22 

INTAKE  AIR  TEMP.  K 

312 

312 

312 

312 

312 

312 

312 

312 

WATER  TEMP.  K 

353 

353 

353 

353 

353 

353 

353 

353 

EXHAUST  TEMP  K 

570 

573 

583 

589 

597 

IE9I 

612 

620 

SMOKE  LEVEL  (Amp.) 

40 

40 

39 

38 

39 

37 

36 

34 

SOOT  AT  STP  (g/m^) 

.698 

.764 

.831 

.935 

.897 

.963 

1.10 

1.35 

SOOT  AT  EXHAUST 
(g/m^) 

.370 

.398 

.428 

.4% 

.452 

.478 

.536 

.716 

NO,  (PPM) 

2300 

2600 

2700 

2000 

2150 

2400 

2500 

NO  (PPM) 

1800 

2000 

2250 

2400 

1750 

2000 

2150 

2300 

FUEL  INJ.  g/min 

10.8 

11.0 

11.4 

11.8 

12.0 

12.0 

12.4 

12.9 

IGNIT.  DELAY  (msec) 

0.64 

0.68 

0.76 

0.60 

0.68 

0.72 

0.76 

INJ.  DURATION  (msec) 

1.24 

1.28 

1.32 

1.36 

1.36 

1.40 

1.48 

1.52 

INJ.  DURATION  (DEG) 

11.2 

11.5 

11.9 

12.2 

12.3 

12.6 

13.3 

13.7 

FUELLING  RATE  *10  ’ 
(g/inj) 

1439 

1465 

1519 

1572 

1598 

1598 

1652 

1718 

IGN.  DELAY  DEG 

5.76 

6.12 

6.84 

7.20 

5.40 

6.12 

6.48 

6.84 

BMEP  (bar) 

3.91 

3.91 

3.91 

3.91 

4.47 

4.47 

4.47 

4.47 

POWER  OUTPUT  (W) 

2199 

2199 

2199 

2199 

2513 

2513 

2513 

2513 

BSFCg/kwh 

294.7 

300.1 

311.0 

321.9 

286.5 

286.5 

2%.0 

308.0 

AIR  FLOW  RATE 
(kg/h) 

21.25 

21.25 

21.25 

21.25 

21.25 

21.25 

21.25 

21.25 

AIR  FUEL  RATIO 

32.80 

32.20 

31.07 

30.02 

29.52 

29.52 

28.56 

27.46 

VOL.  EFFICIENCY  % 

91.38 

91.38 

91.38 

91.38 

91.38 

91.38 

91.38 

91.38 

ROOM  TEMP,  F 

72 

72 

72 

72 

72 

72 

72 

72 
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3.  EXPERIMENTAL  RESULTS 

3.1  PRESENTATION  OF  DATA 

Plots  of  data  for  aU  cases  investigated  are  provided  in  the  appendices.  This 
includes  data  for  fhel  consnmption  (Appoidix  A)»  NOx  (Appendix  B),  soot  (App«dix  C), 
ignition  delay  (Appendix  D)  and  exhaust  gas  temperature  (^pendix  E).  AH  of  the  data 
was  included  in  the  appendices  in  order  to  have  it  organized  and  available  for  scrutiny 
by  interested  parties.  A  page  typteally  contains  four  plots,  r^resenting  four  different 
engine  speeds,  of  the  above  variables  versus  load  or  timing  for  a  common  value  of  tae 
remaining  independent  variable.  Figures  for  hexadecane  fhel  precede  those  for  dodecane 
fuel  in  each  appendix.  The  following  provides  a  discussion  of  general  trmids  observed. 
For  convenience,  certain  example  data  plots  from  the  appendices  are  reproduced  in  the 
following  sections  to  illustrate  these  trends.  Thus  one  can  r^  simply  on  the  examples 
used  here  for  iUustration,  or  refer  to  the  appendices  for  additional  in-depth  evaluation  of 
all  cases. 

3.2  EFFECT  OF  COATINGS  ON  CYLINDER  PRESSURE 

The  effect  of  coating  various  combustion  chamber  surfaces  on  the  cylinder  pressure 
is  illustrated  in  Figure  3.1.  The  four  plots  in  this  figure  are  for  different  injection 
timings.  All  of  the  data  u  for  hexadecane  fiiel,  with  an  engine  speed  of  1500  RPM  and 
a  load  of  447  KPa  BMEP  (the  maximum  load  used  for  these  tests).  In  aU  cases,  the 
baseline  engine  produces  the  highest  cylinder  pressure  followed  in  descending  ordm*  by 
the  coated  head,  all  surfaces  coated,  and  coated  piston  cases.  Some  decrease  in  maximum 
cylinder  pressure  can  be  attributed  to  decreased  volumetric  efficimicy  that  occurs  with 
the  coated  surfaces  (all  tests  were  performed  under  normalfy  aspirated  conditions),  but 
changes  in  the  combustion  process  caused  by  the  surface  coatings  accounts  for  most  of 
the  decrease  as  discussed  in  detail  in  Chapter  4.  This  is  evident  since  operating  the 
engine  with  aD  of  the  surfaces  coated  results  in  the  lowest  volumetric  efficiency,  while 
coating  only  the  piston  crown  consistently  produces  the  lowest  pressure.  Coating  only  the 
piston  effects  the  combustion  process  to  the  greatest  degree  of  all  tiie  insulation  schemes 
and  produces  the  lowest  fuel  consumption  and  NOx  emission  as  described  in  the  following 
sections. 


3.3  EFFECT  OF  COATINGS  ON  BSFC 


Figure  3.2  compares  the  brake  specific  fuel  consumption  (BSFC)  for  the  four 
different  engine  builds  as  a  function  of  load  (BMEP)  for  an  injection  timing  of  20  DBIDC 
for  hexadccane  fticL  Figure  3  J  makes  the  same  comparison  for  dodccanc  focL  Each 
graph  in  these  figures  represents  a  different  engine  speed.  The  injection  timing  valnc  of 
20  DBTDC  was  chosen  arbitrarily  for  purposes  of  these  comparisons.  The  samd‘ trends 
are  present  for  the  other  injection  timings  tested.  The  fhll  set  of  BSFC  data  is  included 
in  Appendix  A.  The  effect  of  varying  the  injection  timing  is  to  increase  BSFC  with 
advancing  timing  as  indicated  in  Figure  3.4.  The  gencrai  trend  exhibited  in  Figures  3.2 
and  3  J  is  that  specific  fuel  consumption  for  the  coated-piston  case  is  better  than  baseline 
for  a  number  of  operating  conditions^  particnlarly  under  low  load  and  at  high  rpm.  The 
coated-head  case  produces  lower  fuel  consumption  for  a  few  combinations  of  operating 
conditions  and  higher  than  baseline  for  others.  When  all  of  the  surfaces  are  coated 
(piston,  head  and  liner),  fhel  consumption  is  always  worse  than  basdine,  probabfy  due  to 
the  reduced  volumetric  efficiency. 

What  appear  to  be  anomalons  readings  in  several  instances,  e.g.  the  low  fuel 
consumption  for  the  coated  piston  case  at  1000  RPM  and  a  load  of  335  KPa  and  the  high 
consumption  for  the  coated  head  case  at  1000  RPM  and  391  KPa,  one  might  at  first 
glance  attribute  to  experimeutal  scatter.  However,  a  review  of  the  entire  set  of  BSFC 
data  in  Appendix  A  show  these  apparent  anomalies  to  be  remarkabfy  consistent  at  the 
same  engine  speed  throughout  the  iujection  timing  range  and  for  both  fheis.  The  strong 
consistency  of  trends  in  the  fuel  cousumptiou  data  adds  coufideuce  to  the  validity  of  the 
measurements.  Uncertainty  in  the  fuel  consumption  data  is  mtimated  to  be  ±10  gr/kwh 
indicated  by  the  uncertainty  bracket  shown  in  Figure  3.2(a).  Apparently,  certain 
combinations  of  operating  conditions  (speed,  timing  and  load)  can  work  together  to 
produce  either  unusually  low  or  uuusually  high  fuel  consumptiou.  The  lowest  overall  fuel 
consumption  for  hexadccane  (270  gr/kwh)  and  dodccane  (260  gr/kwh)  is  produced  by  the 
baseline  en^c  at  2000  rpm,  16  DBTDC  and  a  ioad  of  391  KPa. 
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1000  RPM  &  20  OBTOC 


1500  RPM  &  20  OBTOC 


400 

BMCP  (KPa) 


230 +- 
200 


400 

BMCP  (KPa) 


Fig  3.2(a) 


Fig3J(b) 


2000  RPM  k  20  OBTOC 


2500  RPM  &  20  OBTOC 


250 -t- 
200 


400 

BMCP  (KPa) 


250-1- 

200 


400 

BMCP  (KPa) 


Fig  3  J(c) 


Fig  3.2(d) 


Figure  3.2 


BSFC  Venus  Load  for  Hexadecane  Fuel  at  20  Degrees  BTDC 


BSFC  (gr/kwh)  BSFC  (gr/kwh) 


1000  RPM  &  20  DBTDC  1500  RPM  4  20  D8TDC 


Fig  3.3(a) 


Fig  3.3(b) 


2000  RPM  4  20  OBTOC 


2500  RPM  4  20  OBTOC 


Fig  3.3(c) 


Fig  3.3(d) 


Figure  3.3 

BSFC  Vmus  Load  for  Dodecane  Fud  at  20  Degrees  BTDC 
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Bsrc  ((rn/kwh)  BSTC  (^n/kwh) 


1 000  RPM  k  4.47  Bar 


1500  RPM  k  4.47  Bor  BME 


14  16  18  20  22  24  26 

Timhg  (OBTDC) 


14  16  18  20  22 

Timhg  (DSTDC) 


Fig  3.4(t) 


Fig  3.4(b) 


2000  RPM  k  4.47  Bar  BMEP 


2500  RPM  k  4.47  Bar  BkC 


14  16  18  20  22  24  26 

Timing  (D8T0C) 


14  16  18  20  22  24  26 

Timing  (OBTDC) 


Fig  3.4(c) 


Fig  3.4(d) 


Figure  3.4 


BSFC  Versus  Timing  for  Dodecane  Fuel  at  High  Load 
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3.4  EFFECT  OF  COATINGS  ON  NOx 


Figures  3.5  and  3.6  compare  the  NOx  emissioii  in  PPM  for  each  of  the  four  engine 
bttihis  as  a  function  of  load  for  hcxadccane  and  dodccane,  respectiveiy.  Again,  there  are 
separate  plots  for  each  engine  speed,  and  the  example  data  shown  is  for  a  timing  of  20 
DBTDC.  Complete  data  plots  for  all  runs  are  provided  in  Appendix  B.  UnMitainty  in 

the  NOx  data  is  estimated  at  ±50  PPM,  approximately  represented  by  the  size  of  the  data 
symbob  in  the  figures. 

Ceramic  coating  the  pbton  crown  has  a  very  favorabb  effect  on  the  NOx  emission 
of  thb  engine.  The  decrease  in  NOx  for  hexadecane  at  20  DBTDC  timing  ranges  from 
zero  percent  at  high  speed  and  low  load  to  25  percent  at  low  speed  and  high  load.  The 
decrease  under  high  load  conditions  ranges  from  20  to  34  percent  depending  on  injection 
timing,  regardless  of  engine  speed.  Possibb  reasons  for  the  decreased  NOx  production 
are  rebted  to  the  effect  that  the  hotter  pbton  surface  has  on  the  combustion  process  u 
discussed  in  the  next  chapter  in  conjunction  with  the  KIVA  modeling  results.  Lowered 
NOx  emission  coupled  with  improved  fuel  economy  would  appear  to  make  coating  of  the 
pbton  crown  a  desirable  engine  modification. 

Coating  of  the  cylinder  head  by  itself,  on  the  other  hand,  does  not  significantly 
reduce  NOx  emission  and  under  some  operating  conditions  causes  an  increase.  Thb  b 
especially  true  with  dodecane  fuel,  where  coating  the  head  causes  increased  NOx  under 
most  conditions.  The  unimproved  or  degraded  NOx  emission,  coupled  with  simUarfy 
unimproved  or  degraded  fuel  consumption,  provide  little  incentive  to  ceramb  coat  the 
cylinder  head. 

Operating  the  engine  with  coatings  on  all  combustion  chamber  surfaces,  Le.  pbton, 
head  and  cylinder  liner,  produces  NOx  emission  values  that  fall  between  those  for  the 
pbton  alone  and  for  the  head  alon^Xx  One  might  expect  thb  result  based  on  the 
reasoning  that  the  benefidal  effects  of  coating  the  pbton  are  ofbet  somewhat  by  the  often 
adverse  effects  of  coating  the  head.  In  addition,  the  insnbting  effect  of  coating  all  three 
surfaces  results  in  slightly  higher  cylinder  gas  temperatures  throu^out  the  combustion 
cycle  which  will  tend  to  increase  NOx  as  shown  by  the  KIVA  modeling  resulb  discussed 
in  Chapter  4.  Figure  3.7  shows  the  effect  of  varying  injection  timing  on  NOx  emission 
at  high  load  and  four  engine  speeds  for  dodecane  fueL  The  trend  b  for  ihereasing  NOx 
with  advancing  injection  timing  as  would  be  expected  due  to  the  increased  ignition  delay 
and  premixed  combustion  fraction  resulting  in  higher  peak  cylinder  pressure  and 
temperature. 
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Figure  3.5 

NOx  Versus  Load  for  Hexadecane  Fuel  at  20  Degrees  BTDC 
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Figure  3.7 

NOx  Versus  Injection  Timing  for  Dodecane  Fuel  at  High  Load 
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3.5  EFFECT  OF  COATINGS  ON  IGNITION  DELAY  PERIOD 


Figures  3.8  and  3.9  show  examples  of  ignition  delay  versus  injection  tuning  at  light 
load  for  hexadecane  and  dodecane  fuel  respectively.  Data  for  all  cases  b  provided  in 
Appendix  C.  Ignition  delay  was  measured  with  the  digital  waveform  analyzer  used  to 
collect  the  data  by  placing  the  cursor  at  the  point  on  the  needle  lift  trace  where  the  needle 
opens  and  then  advancing  to  the  point  on  the  pressure  trace  where  rapid  pressure  rise 
begins.  At  the  maximum  25  KHz  data  rate  of  the  instrument,  the  time  between  data 
points  b  0.04  ms,  and  thb  b  assumed  to  represent  the  uncertainty  in  the  ignition  delay 
measurements.  The  uncertainty  of  ±0.04  ms  b  indicated  by  the  brackets  in  the  figures. 

The  ignition  delay  values  for  both  fueb  lie  between  0.4  and  1.1  ms  depending 
mainly  on  speed  and  injection  timing,  and  to  a  lesser  extent  on  load.  These  values  are 
typical  for  high  cetane  fiieb. 

At  low  speeds  the  coated  pbton  appears  to  decrease  the  ignition  deby  relative  to 
the  baseline  engine  as  predicted  by  the  KIVA  resulb  discussed  in  the  next  chapter.  The 
shortened  ignition  delay  b  the  reason  proposed  in  Chapter  4,  based  on  computer 
modeling,  for  reduced  cylinder  gas  temperature  and  pressure  and,  consequently,  reduced 
NOx  relative  to  the  baseline.  In  many  cases  the  values  of  ignition  delay  for  the  coated- 
pbton  case  and  the  baseline  case  are  within  the  limib  of  uncertainty  and  are  therefore 
inconclusive.  In  a  few  cases,  the  measured  ignition  delay  for  the  coated-pbton  case  b 
clearly  longer  than  that  for  the  baseline  case,  e.g.  for  hexadecane  at  2500  RPM  and  335 
KPa  BMEP  (Figure  C2(d),  Appendix  C).  For  thb  same  case,  the  NOx  emission  of  the 
coated-pbton  engine  b  slightly  less  than  the  baseline  engine,  in  contradiction  to  our  theory 
for  decreased  NOx.  The  conclusion  that  can  be  drawn  from  thb  b  that  ignition  delay  b 
not  the  only  controlling  parameter  in  the  NOx  reaction  and  additional  combustion 
phenomena,  of  which  we  do  not  have  a  complete  understanding  at  thb  time,  are  affected 
by  the  ceramic  coatings.  The  lack  of  in-cylinder  gas  temperature  measurements  hinders 
our  ability  to  completely  understand  the  effect  of  coatings  on  the  total  combustion 
process. 

For  the  cases  of  coated  head  alone  and  all  surfaces  coated,  no  unequivocal 
statements  can  be  made  concerning  the  effect  on  ignition  delay.  For  hexadecane  fuel, 
coating  the  head  produces  the  longest  ignition  delay  values  under  low  load  conditions, 
while  coating  all  the  surfaces  produces  the  longest  values  under  high  load  conditions.  For 
dodecane  fuel,  coating  all  surfaces  generally  produces  the  longest  ignition  delays 
regardless  of  load. 
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Figure  3.8 


Ignition  Deiay  Versus  Timing  for  Hexadecane  Fuel  at  223  KPa  BMEP 
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Figure  3.9 


Ignition  Delay  Versus  Timing  for  Dodecane  Fuel  at  223  KPa  BMEP 


3.6  EFFECT  OF  COATINGS  ON  SOOT  EMISSION 


Soot  concentration  values  in  gr/m’  were  caknlated  from  the  photocell 
output  of  the  USPHS  smoke  opacity  meter  by  calibrating  the  meter  output  with 
a  set  of  neutral  density  filters  and  using  the  concentration  versus  opacity 
data  from  Reference  11.  The  measurement  uncertainty  for  the  soot 
concentration  values  is  unknown,  but  the  data  is  useful  for  observing 
relative  effects  of  the  different  coating  schemes  on  soot  production. 

Figure  3.10  u  representative  of  most  of  the  soot  data.  It  provides  soot 
concentration  versus  load  for  heudecane  fuel  at  22  DBTDC  injection  timing 
for  four  engine  speeds.  Most  of  the  trends  displayed  in  this  figure  hold  true 
for  the  rest  of  the  soot  data  which  are  presented  in  Appendix  D.  Coating  all 
of  the  combustion  chamber  surfaces  has  a  deleterious  effect  on  soot  emisuon, 
particularly  at  lower  speeds  and  higher  loads.  The  effect  on  soot  emission 
of  coating  the  piston  alone  depends  on  engine  speed.  Coating  the  piston 

results  in  higher  soot  production  tlmn  baseline  for  some  engine  speeds  and 

lower  at  others.  Likewise,  the  effect  of  coating  the  head  alone  depends  on 
engine  speed.  At  1500  RPM  the  coated  head  produces  relatively  low  soot 
emissions,  while  at  2500  RPM  it  produces  relatively  high  amounts  of  soot. 

Possible  explanations  of  these  results  based  on  computer  simulations 
are  not  available  at  this  time.  The  soot  modeling  described  in  the  following 
chapter  has  so  far  been  applied  only  to  the  baseline  engine.  While  the  results 
for  the  baseline  case  are  encoun^ing,  the  soot  model  has  not  yet  been 

implemented  in  conjunction  with  the  coated  surface  temperature  model  in  order 
to  predict  the  effects  of  various  surface  coatings  on  soot  emission.  This 

work  is  continuing  and  will  be  published  in  technical  paper  form  when 
completed. 
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Figure  3.10 


Soot  Concentration  Versus  Load  for  Hexadecane  at  22  DBTDC 


3.7  EFFECT  OF  COATINGS  ON  EXHAUST  GAS  TEMPERATURE 


Typical  exhaiut  gas  temperature  data  are  shown  in  Figure  3.11  which  shows 
exhaust  temperature  versus  load  for  hexadecane  at  20  DBTDC  injection  timing.  The 
same  trends  are  exhibited  by  the  rest  of  the  exhaust  temperature  data  provided  in 
Appendix  E.  From  these  figures  several  interesting  phenomena  can  be  noted.  Coating 
the  head  by  itself  has  little  effect  on  the  exhaust  temperature  relative  to  basdhie  except 
at  low  speed  where  it  actually  reduces  exhaust  temperature.  This  trmid  is  predicted  by 
the  KIVA  results  (e.g.  see  Figure  4.27).  A  study  by  Morel,  et  aL  [12]  indicated  that  for 
a  non-insulated  diesel  engine,  49  percent  of  the  heat  loss  from  the  hot  cylinder  gas  occurs 
through  the  piston,  32  percent  occurs  through  the  head  and  the  remaining  19  percent 
occurs  through  the  liner.  Thus,  insulating  the  head  should  have  less  impact  on  the  total 
heat  transfer  and  the  exhaust  gas  temperature.  Peak  cylinder  gas  temperatures  are  less 
than  baseline  due  to  the  effect  of  the  hi^er  surface  temperature  on  the  combustion 
process,  as  discussed  in  Chapter  4. 

The  exhaust  temperature  for  the  coated-piston  case  is  typically  higher  than  baseline 
due  to  the  reduced  heat  transfer  through  the  piston  to  the  oiL  Since  approximately  half 
of  the  heat  transfer  takes  place  through  the  piston,  adding  a  low  conductivity  coating  to 
the  piston  crown  should  significantly  reduce  heat  transfer  and  increase  exhaust 
temperature. 

Coating  all  of  the  surfaces  produces  the  highest  exhaust  temperatures  as  would  be 
expected  due  to  the  larger  reduction  in  heat  loss  during  the  expansion  stroke  when  all 
surfaces  are  insulated. 
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Figure  3.11 

Exhaust  Gas  Temperature  Versus  Load  for  Hexadecane  at  20  DBTDC 
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4.  COMPUTER  MODELING 


4.1  BACKGROUND 

The  KIVA-11  code,  developed  et  Loo  AluBoe  NetioBel  Laboratory  [13],  toKee 
dther  the  2>D  or  the  3-D  oasteady  eqaatioBt  of  aiotioB  of  a  chMBkaby  reactive  adxtara 
of  ideal  gaaes  incladiag  the  dyaaadcs  of  a  Uqaid  fhel  spnqr  aad  the  cooplii^{  betweea  the 
spray  aad  the  gas.  The  code  was  developed  with  appheatioB  to  iateraal  coaibartioa 
eagfaies  specifically  ia  miad,  aad  coataias  a  aaaibcr  of  featares  to  CMitttate  this 
appUcatioB,  sach  as  gas  flow,  Uqaid  faei  iajectioa,  spray  dyaaasics,  evaporatioa,  heat 
traasfer,  combastioa,  species  traasport,  aad  aifada^ 

KIVA-n  was  chosea  for  this  work  becaase  of  Ae  docaaieatatioB  Aat  is  available 
[13,14,13]  aad  the  fact  that  it  is  bdag  widely  ased  by  oAer  researchers  ia  iadastry, 
acadeada  aad  natioBal  laboratories  to  modd  combastioa  aad  easissioas  ia  boA  sparii 
igaitioa  aad  diesel  eagiaes,  e.g.  Aaisdea  et  aL[16],  0*Roark  aad  Amsdea  [17],  Gmitry  et 
aL  [18],  Naber  aad  Reitz  [19],  Kao  aad  Rdtz  [20],  Reitz  and  Dhraker  [21],  PAchon  [22], 
Zellat  [23],  Gibson  et  aL  [24],  Varaavas  aad  Assaais  [25].  Ia  additioa,  Ramos  [26],  aad 
Marfcatos  [27]  docameated  KIVA  A  their  recent  books. 

The  origAal  KTVA-II  code  emptoys  the  ArrfaraAs  sA^  reactAn  combastAa 
model  (SR  model)*  This  modd  is  saiAbk  for  premAed  combastAa  systems  sach  as  spark 
igaitAa  eagAcs.  However,  A  dksd  eagAes  the  combastAa  process  is  maialy  difAsAaal 
after  a  short  preaiAed  period.  The  mixiag-coatroUed  eddy-break-op  modd  (EBU  model) 
devdoped  by  Magaassea  et  aL  [28-31]  provides  a  better  physical  reprcscaAtioa  of  Ac 
diffosAaal  combastAa  process. 

A  this  sAdy,  lOVA-H  was  modified  to  Aclade  Ac  EBU  combastAa  model,  aad 
resalts  from  boA  modeb  are  compared  with  experiments.  New  model  constants  were 
Atrodaced  fbr  Ac  origAal  Arrhenias  combastioa  modeL  These  constants  were  ased  for 
different  combinations  of  Aad,  AjectAn  timAg  and  speed  and  were  foand  to  give  good 
agreement  between  the  predicted  combastAa  and  experimento  for  both  hexadccane  and 
dodccanc  AcL 

As  part  of  this  work,  a  soot  model  was  also  added  to  KIVA-IL  The  model  for  the 
rate  of  soot  formatAn  developed  by  Tesner  et  aL  [3233]  was  implemented  to  calcnAte  the 
amonnt  of  soot  formatAn,  and  Ac  model  of  soot  combnstion  developed  by  Magnussen 
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asd  mcrtagur  |28|  coaccraiag  the  behavior  of  soot  ia  torbakat  flaaict,  was  iasidcamted 
to  coaipatc  Ac  oaddatioa  rate  of  soot 


4.2  GENERAL  STRUCTURE  OF  KIVA-n  CODE 

The  ga*>phase  sohitioa  procedare  ia  iOVA-D  is  based  oa  a  fiaite  vohuae  method 
called  the  ALE  (Arbitraiy  Lagraagiaa-Ealeriaa)  method  which  fuilitates  cakalatioos 
with  carved  aad  chaagiag  booadaries.  Spatial  differeaces  arc  formed  oa  a  flaitc- 
differcacc  mesh  Aat  sabdividcs  the  compatathaud  regioa  iato  a  aamber  of  saudl  cells  that 
are  hexahedroas.  The  code  fcatorcs  a  Stochastic  Particle  TcchaiquCf  aa  cilieieat  aad 
accurate  method  for  solviag  the  spray  dyaaoilcs,  based  oa  the  Monte  Cario  aad  discrete 
particle  methods.  The  traasicat  sohitioa  is  aiardied  oat  ia  a  scqueacc  of  flaitc  time 
iacremcats  called  cycles  or  timesteps.  Oa  each  cycle  Ac  values  of  Ac  dcpcadoit  variables 
are  cakulated. 

The  KIVA'II  computer  program  coasists  of  a  set  of  sabroatiacs  coatroUed  by  a 
short  mala  program.  It  was  wiittea  specifically  for  use  oa  the  CRI  Cray  fiuaily  of 
computers,  operating  uader  Ac  Cray  Time  Sharing  System  (CTSS)  and  odng  Ae  Cray 
FORTRAN  (CFT  and  CFT77)  compilers. 

AU  compuAtional  runs  were  made  on  the  CRAY  Y-MP  at  Ae  NorA  Carolina 
Supercomputiag  Center  (NCSC).  KTYA-D  contains  statements  peculiar  to  Ac  CFT  and 
CFT77  compilers  that  permit  vcctorization  of  many  of  Ac  loops  ia  the  sabroatiacs.  The 
vectorization  of  Ac  code  improved  Ac  ran  time  by  a  factor  of  five.  A  typical  ran  for  Ac 
Ricardo  Hydra  engine  nmoiation,  including  the  soot  model,  begiaiag  at  138  degrees 
BTDC  and  ending  at  510  degrees  ATDC,  required  about  1200  CPU  seconds  oa  Ae  Cray 
Y-MP. 


4.3  COMPUTATIONAL  PARAMETERS 

For  simplicity,  a  2-D  axisymmctric  simulation  of  Ae  combustion  in  Ac  Ricardo 
Hydra  direct  injection  diesel  engine  was  used.  Figure  4.1  shows  Ac  KIVA-Il 
compuAtional  mesh  of  Ac  Ricardo  Hydra  combustion  chamber  at  138**  BTDC  (point  of 
intake  valve  closure),  and  again  at  TDC.  The  mesh  has  23  cells  in  Ac  radial  direction,  1 
cell  in  Ac  azimuAal  dircetton  and  24  cells  in  the  axial  direction  (nx»23,  ny^l,  iiz»24). 
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At  TDC 


Figure  4.1 

Computational  Mesh  Used  to  Model  Ricardo  Hydra  DI  Engine 
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The  calculations  were  started  when  the  intake  valve  ch>sed  and  ended  at  90  degrees 
crank  angle  ATDC.  The  initial  species  was  air  which  was  assumed  to  he  21  percent 
oxygen  and  79  percent  nitrogen.  The  initial  densities  of  the  oxygen  and  nitrogen  were 
assigned  so  that  the  initial  pressure  in  the  combustion  chamber  as  the  intake  valve  closed 
was  96  Kpa.  The  initial  air  temperature  was  assumed  to  be  310  K. 

To  calculate  the  dynamics  of  the  spray,  the  model  requires  the  radii  and  velocities 
of  the  droplets  at  the  injector  (15].  The  magnitude  of  the  injection  velocity,  Vj,j,  is 
determined  from 


p,it 


(4.1) 


where: 


is  the  injected  fuel  mass,  g/cycle. 
is  the  fuel  density,  g/cm^ 
is  the  injector  noazle  radius,  cm 
is  the  injection  duration,  s. 

The  sauter  mean  radius  was  assumed  to  be  9.5  micron.  The  model  requires  the 
injection  angle,  injected  fuel  mass,  fuel  density,  and  injection  duration.  These  values  were 
taken  from  the  experimental  data. 
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4.4  COMBUSTION  MODELS 


4.4.1  Single  Reaction  Model 

The  single  reaction  Arrhenius  modd  is  the  standard  conbostion  model  built  into 
the  KIVA-n  code.  Oxidation  of  the  fhel  b  assumed  to  be  modeled  by  a  one-step  global 
kinetics  scheme.  The  complete  set  of  chemical  reactions  b  given  below. 

Two  different  fiieb  were  used  in  both  the  experiments  and  the  modeling,  vb 
hexadecane  and  dodecane  The  oxidation  reaction  for  the  two  fiicb  are 

given  by: 

*/i 

2  +  49  O,  32  CO,  +  34  ff,0 


and 


2  +  37  24  CO,  + 


26  H^O 


(4J) 


The  reaction  rate  for  these  oxidation  reactions  b  given  by: 

(0  -  [Fuel\*  {Oxidizer]^ 

The  coefficient,  ,  b  assumed  to  be  of  a  generalized  Arrhenius  form, 

h  -  4.  ^  “P 


(4.4) 


(4.5) 


where: 

Af  b  the  Arrhenins  coefficient 
^  b  the  forward  temperature  exponent 
Ef  b  the  activation  temperature,  K. 

The  most  common  choice  of  parameters  in  the  global  rate  expression  b  a  =  1  and 
b  -  1,  aMnming  that  the  overall  reaction  b  approximately  first  order  with  respect  to  both 
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fuel  and  oxidizer.  The  other  parameters,  /i,  and  E,  ,  can  be  adjusted  so  that  the 

predicted  Pressure-Crank  Angle  diagram  matches  the  experimental  value.  It  was  found 
in  this  study  that  the  best  agreement  with  experiments  was  obtained  with  values  of  the 
fuel  exponent  a  -  0.25  and  the  oxidizer  exponent  b  -  1.5  for  both  hexadecane  and 
dodecane.  The  best  values  of  the  Arrhenius  coefficient  and  activation  temperature  that 

gave  satisfactory  results  for  both  fuels  were  found  to  be  A,-9.0eid  and  E,~l.Se4  K 

respectively. 

4.4.2  Eddy-Break-Up  Model 


In  diffusion  flames,  fuel  and  oxygen  occur  in  separate  eddies.  Because  the  chemical 
reactions  in  most  cases  are  very  fast,  it  can  be  assumed  that  the  rate  of  combustion  will 
be  determined  by  the  rate  of  intermixing  of  fuel  and  oxygen  eddies  on  a  molecular  scale: 
in  other  words,  by  the  rate  of  dissipation  of  the  eddies.  Because  fiiel  and  oxygen  appear 
as  fluctuating  intermittent  quantities,  there  will  be  a  relationship  between  the  fluctuations 
and  the  mean  concentration  of  the  species.  Consequently,  the  rate  of  dissipation  can  be 
expressed  by  the  mean  concentration  of  the  reacting  species.  The  EBU  model  gives  a 
physically  based  representation  of  combustion,  because  it  takes  into  account  the  effect  of 
turbulence  on  the  mean  chemical  reaction  ratesHx  In  this  study,  the  model  of  Magnussen 
and  Hjertager  was  used  [28].  This  model  relates  the  rate  of  combustion  to  the  rate  of 
dissipation  of  turbulent  eddies  and  expresses  the  rate  of  reaction  by  the  mean 
concentration  of  a  reacting  species,  the  turbulent  kinetic  energy,  and  the  rate  of 
dissipation  of  this  energy.  This  model  differs  from  the  EBU  model  of  Spalding  [3435]  in 
relating  the  dissipation  of  eddies  to  the  mean  concentration  of  intermittent  quantities 
instead  of  the  concentration  fluctuations.  This  is  advantageous  in  view  of  the  uncertainty 
in  determining  the  concentration  fluctuations  of  the  reacting  species. 

Accordingly,  in  the  case  where  the  adequate  presence  of  fuel  is  the  controlling 
factor,  the  rate  of  combustion  of  fuel  can  be  expressed  by: 


A  ,  Cj . 


(4.6) 


where: 

A  is  a  constant 

Cf  is  the  local  time-mean  fuel  concentration 

k  is  the  turbulent  kinetic  energy  cm^/s^. 

8  is  the  rate  of  dissipation  of  turbulent  kinetic  energy  emVs^ 


32 


In  regions  of  the  flame  where  the  time  mean  concentration  of  fuel  is  h^  and  the 
oxygen  concentration  b  low,  oxygen  will  be  the  reacting  species  that  controb  the  rate  of 
combustion.  Accordingly,  the  rate  of  combustion  in  thb  case  can  be  expressed  by: 


(4.7) 


where: 

is  the  local  time-mean  oxygen  concentration 
r,  b  the  stoichiometric  oxygen  requirement  to  bum  1  kg  fiieL 

In  premixed  turbulent  flames,  fuel  and  oxygen  will  occur  in  the  same  eddies.  These 
eddies  will  be  separated  by  eddies  containing  hot  combustion  products.  The  rate  of 
combustion  will,  in  thb  case,  be  determined  by  the  same  mechanbm  described  above. 
However,  an  extra  equation  which  takes  care  of  the  dissipation  of  the  hot  eddies  most  be 
added  in  cases  where  the  concentration  of  hot  combustion  product  b  low.  The  combustion 
rate  can  then  be  expressed  by: 

k.-A-B-Si—.  [1  1  (4^ 

'  !./>  [*j 

where: 

B  is  a  constant 

Cp  b  the  local  time-mean  concentration  of  reaction  products. 


Thus,  the  EBU  model  b  applicable  to  diffusion  as  well  as  to  premixed  flames.  The 
equation  that  gives  the  lowest  reaction  rate  b  the  one  that  determines  the  local  rate  of 
combustion.  Therefore,  combining  equations  (4.6-4.8),  the  rate  of  combustion  can  be 
expressed  by: 


-  A  .  min 


P 

</  *  *  *> 


8 

1 


(4.9) 


For  the  model  constants,  Magnussen  suggested  the  values  of  ^4  -4.0  and  B  =  0.5, 
but  Gosman  [36]  used  A  =20  and  B  =  2.5,  while  Pinchon  [22]  used  A  =  16  and  B  =  2  and 
Varaavas  [25]  used  A  =  0.5  and  B  =  0.5.  It  was  found  in  thb  study  that  the  original  values 
used  by  Magnussen  (A  =  4,  B  0.5)  gave  the  best  overall  agreement  with  the 
experimental  data. 
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The  EBU  model  has  no  provision  for  initiating  combustion.  To  circumvent  this 
problem,  the  SR  chemical  kinetic  model  was  employed  up  to  the  point  of  auto-ignition 
after  which  the  calculation  procedure  was  switched  to  the  EBU  model  Switching  from 
the  chemical  kinetics  model  to  the  EBU  model  depends  on  whether  the  ratio  of  the  local 
turbulent  time  scale,  *  to  the  local  kinetic  time  scale,  r, ,  called  the  Damkohler  number 
[23]  is  less  than  or  greater  than  unity. 

When  T^l  T,<  1,  the  process  is  controlled  by  chemistry  since  the  chemical  reaction 
time  scale  is  much  larger  than  the  turbulent  time  scale. 

When  T„  I  T,  >  1,  the  process  is  mixing  controlled  since  the  time  required  for 
chemical  reactions  is  negligible  compared  to  the  time  required  for  mixing 

The  two  time  scales  are  defined  by  the  relations: 


and 


T,  -  ^  .  p  txp(-E/RJ) 


where: 

is  the  turbulence  mixing  time  scale 
X,  is  the  chemical  reaction  time  scale 
A  is  the  pre-exponential  coefficient  in  the  kinetics  formula 
p  is  the  density 

E  is  the  activation  energy 

R,  is  the  universal  gas  constant 
k  is  the  turbulent  kinetic  energy 

8  is  the  dissipation  rate. 

It  was  decided  to  switch  from  the  SR  to  the  EBU  model  of  combustion  after  the 
ignition  delay  period  when  the  combustion  becomes  controlled  by  mixing  The  cylinder 
gas  temperature  was  used  to  trigger  the  switch.  A  cell  gas  temperature  of  1500  K  was 
choosen  to  signal  the  switch  to  the  EBU  model.  This  cell  temperature  is  generally 
exceeded  within  one  time-step  of  the  onset  of  combustion.  By  contrast,  the  maximum 
motoring  temperature  (no  combustion)  is  about  950  K. 
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Although  there  is  much  evideuce  to  support  the  dependmce  of  the  combustwn  rate 
on  turbulent  mixiug,  it  is  important  to  note  that  the  mixing  controlled  EBU  model  is 
inadequate  to  predict  combustion  near  walls,  self-ignition,  pollutant  formation,  and  lean 
and  rich  flammability  limits.  These  phmiomena  are  controlled  by  chemical  kinetics. 

4.4.3  Results  of  Combustion  Modeling 

Comparisons  of  the  P-6  diagrams  obtained  for  the  baseline  engine  with  both  the 
Eddy-Break-Up  model  and  the  single  reaction  model  with  experimental  results  have  been 
made  for  several  different  loads  and  injection  timings,  and  some  representative  results  are 
presented  here.  As  seen  in  Figures  4.2  and  4J  for  hexadecane  fuel  under  two  different 
loads  for  the  same  timing  and  engine  speed,  the  cylinder  pressure  rises  faster  (rate  of  heat 
release  is  greater)  during  premixed  combustion  for  the  SR  Model  compared  to  the  EBU 
ModeL  Beginning  at  the  point  of  ignition,  the  rate  of  heat  release  shown  by  the  EBU 
Model  is  representative  of  diffusion  combustion.  There  appears  to  be  no  premixed 
pressure  spike.  The  EBU  Model  predicts  slightly  lower  peak  pressure  than  the  SR  Model 
in  all  cases  examined.  A  similar  trend  b  noted  for  dodecane  fuel  as  well,  as  shown  in 
Figures  4.4  and  4.5. 

It  b  abo  evident  from  the  P-6  diagrams  that  the  predicted  ignition  delay  b  always 
longer  than  the  experimental  value  by  approximately  3-5  degrees.  Recall  that  the  EBU 
Model  has  no  mechanbm  for  autoignition,  and  that  when  using  the  EBU  Model  we  rely 
on  the  Arrhenius  SR  Model  to  initiate  combustion.  We  thus  see  identical  ignition  delay 
periods  for  both  modeb.  Ignition  delay  in  the  SR  model  b  primarily  a  function  of  the 
activation  temperature  constant,  E,,.  The  choice  of  activation  temperature  that  gave  the 
best  overall  results  for  the  SR  Model  over  a  wide  range  of  operating  conditions,  i.e. 
En=1.5*10\  resulted  in  slightly  longer  ignition  delay  periods. 

Figures  4.6  and  4.7  show  the  predicted  cylinder  pressure  traces  for  three  different 
injection  timings  for  the  EBU  Model  and  SR  Model,  respectively.  The  maximum  cylinder 
pressure  increases  with  advancing  injection  timing  as  it  should,  due  to  a  larger  proportion 
of  premixed  combustion. 

Figure  4.8  shows  the  predicted  cylinder  gas  temperature  versus  crank  angle  for 
both  combustion  modeb.  The  EBU  Model  predicts  slower  gas  temperature  rise  and  a 
lower  peak  gas  temperature  than  the  SR  Model,  in  keeping  with  the  lower  peak  cylinder 
pressures. 
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CYUNOER  PRESSURE  (ATM) 


CRANK  AHGL£  (0€C)  CRANK  ANGLE  (DEO) 


Flf«K  4.6 

CyiMcr  Pranrc  Vcmu  Craak  Aagle  Far  Hciadacaac  Fatl 
EBU  Modd  for  DMIitrcal  lajadioa  Tiailag* 


FIfarc  4.7 

CjrtMcr  Proaorc  Vcfaaa  Craak  Aagle  tor  Hcaadccaac  Fad 
SR  Modd  tor  DUtorcat  laJccMoa  Tlaiiagi 


Crank  Anglo  (DEC) 


Figar*  4J 

CyMadcr  Tcapcratarc  Verm  Craak  Aagle  For  Hcaadccaac 
Slagle  Rcacifoa  Modd  Verm  Eddy>Brcak*Dp  Modd 


37 


4.5  NTTRIC  OXIDE  MODELING 


4.5.1  NO  Model 

Nitric  oiidc  (NO)  is  formed  in  diesel  engines  by  the  high  temperature  oxidation  of 
atmospheric  nitrogen.  In  oxygen-rich  (lean)  regions,  die  elementary  process  for  the 
formation  reaction  of  NO  b  expressed  as: 

2N^  +  ^  2N  f  2NO 

20,  +  N,  ^  20  ^  2NO 

In  excess-fuel  (rich)  regions,  the  oxygen  concentration  required  in  equation  (4.10) 
b  low,  and  the  reaction  of  N  and  OH  (which  b  formed  during  the  decomposition  of  the 
fuel)  becomes  important  Thb  reaction  b  given  by: 

N,  +  20H  ^  2NO  +  IH  (^12) 


(4.1f) 

(4.11) 


The  elementary  reaction  processes  given  by  equations  (4.10)  and  (4.11)  are  known 
as  the  Zeldovich  mechanbm  When  the  elementary  reaction  process  represented  by 
equation  (4.12)  b  further  add^,  the  total  reaction  b  called  the  extended  Zeldovich 
mechanbm  [38]. 


Assuming  a  steady  state  approximation  for  the  nitrogen  atom  concentration,  the 
rate  expression  fits  the  extended  zeldovich  mechanbm: 


d  [M>] 
dt 


-  2  {0\  [N,] 


1  -  [NO]*  /  K  [O,]  [iV,I 
1  +  [NO]  /  (  [0,1  +  [OH]  ) 


(4.13) 


where: 

[  ]  b  tte  species  concentration 
k.  b  the  forward  reaction  coefficient  for  the  reaction  r 

Jr 

k.  b  the  backward  reaction  coefficient  for  the  reaction  r 

•r 

K  b  the  dissociation  equilibrium  constant  =  i  )(  ^/,  ^  )  * 


The  reaction  rate  constants  used  in  KIVA-Il  are  [39]: 
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A,  -  15587*10*^  exp(-67«7/7) 
ik^  -  13  •10“ 

-  2.6484*10'*  exp( -59418/7) 

-  1.6*10'*  cxp(19678/7) 
k,  -  2.123  *10'*  exp( -57020/7) 

A 

A.  -  0.0 

•4 

The  required  concentrations  were  computed  from  die  foUowinf  set  of 
equilibrium  reactions: 

^  2H 
O,  20 
^  2N 
^  20H 

O,  +  2Hfi  ^  40H 
2  CO  +  O,  ^  2CO, 

The  following  relation  was  used  to  convert  the  NO  values  in  gm/cycle,  as 
calculated  by  KlYA,  to  concentration  in  particles  per  million  (ppm)  for  comparison  with 
experiments, 
where: 

is  the  total  mass  flow  rate  (air  +  fhel  injected) 
rps  is  the  revolution  per  second. 
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4.5.2  ResoHs  of  NO  Modeliog 


Figures  4.9  and  4.10  coaspare  the  NO  crnkyon  versus  crank  angle  for  botii  nsodds 
with  the  expcriniMital  values  for  hexadccanc  ftiel  and  two  different  injection  timings.  The 
EBU  model  predkts  lower  NO  emissions  than  the  SR  model  due  to  the  tower  peak 
cylinder  temperature.  The  values  predicted  by  the  SR  model  are  in  hotter  agrecmMt 
with  the  experimental  values,  although  they  arc  stiU  somewhat  below  the  measured  values. 
Figures  4.11  and  4.12  show  similar  NO  emission  results  for  dodecanc  hicL  For  dodccanc, 
the  SR  model  gives  very  good  agreement  with  the  measured  NO  values. 

Figures  4.13  -  4.15  show  the  effect  of  injection  timing  on  NO  emission  for 
hexadccanc  for  three  different  toads,  comparing  experimental  results  with  both  models, 
th  models  nnderpredict  NO,  but  produce  the  correct  trend  of  increasing  NO  emtoston 
ih  advanced  timing.  Figure  4.16  shows  similar  data  for  dodecanc  at  tow  toad. 
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NO  Emisstons  (PPM) 


FIgara  4.9 

NO  EaiHtM  Venn  Craak  Aagli  For  H«ui4acaac  Fad 
ISM  RFM,  1«  Deg  VrOC,  m  KPa  BMET 


PIgan  4.19 

NO  EMMaa  Venal  Cr^  AagM  tar  Htudccaac  Fad 
ISM  RFM,  M  Di«  nVC,  Ml  KFa  HNEP 


FIgarv  4.12 

NO  Enlnioa  Vcnai  Craak  Aagk  tor  Dadceaac  Fad 
ISM  RPM,  U  D(«  VIDC,  tt3  KPa  BMEP 


FIgara  4.11 

NO  Eaiinioa  Venw  Craak  Aagk  Far  Dadceaac  Fad 
ISM  RPM,  Id  Deg  nDC,  223  KPa  BMEP 
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(PPM)  NO  EMSSUNS  (PPM) 


MJCCnON  1MM40  (OCO-BTOC) 


Figarc  4.13 

NO  EmMm  VcrfM  laJcrtiM  TMa|  For  Hcudccaac  Fad 
ISM  RFM.  233  KPa  IMEP 


Fipir*  4.14 

NO  EMMm  Vtnw  mtdlM  1taii«  Far 
ISM  RFM,  338  KPa  MIIP 


Faai 


INJECTION  THJMO  (OEG-BTOC)  INJECTION  TMNO  (DEO-BTOC) 


Figare  4.15 

NO  EnMoa  Vcnaa  laJcctloB  Tiaiiag  For  Hraadccaac  Fad 
ISM  RFM,  391  KFa  BMEF 


Flgan  4.16 

NO  EnMoa  Vcraot  la|cctioa  Tiariag  for  Dodccaoc  Fad 
ISM  RFM,  223  KFa  BMEF 
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4.6  SOOT  MODELING 


4.6.1  Review  of  Existing  Soot  Models 

Most  modeb  appearing  in  the  literature  for  the  formation  and  oxidation  of  soot 
are  empirical  and  much  of  the  data  used  in  the  development  of  the  modeb  are  based  on 
laboratory  burner  flames  rather  than  diesel  engines.  The  literature  on  soot  has  been 
reviewed  by  several  researchers,  including  Haynes  and  Wagner  [40],  and  most  recently 
by  Morel  and  Keibar  [41].  Although  many  useful  observatwns  have  been  made,  no 
generally  applicabb  soot  emission  model  has  yet  been  devbed. 

Khan,  et  aL  [42,43]  presented  a  soot  formation  correbtion  based  on  high  pressure 
diesel  engine  data.  The  parameters  used  in  thb  model  are  local  values  of  temperature, 
unbumed  hydrocarbon  concentration,  and  equivalence  ratio.  A  corresponding  soot 
combustion  scheme  was  not  suggested.  The  rate  of  soot  formation  was  described  as: 

—  s  exp(-2000/r)  (4.15) 

dt  ’ 


where: 

s  b  the  soot  mass 

b  the  equivalence  ratio  in  the  formation  zone 
Pf  b  the  partial  pressure  of  fuel  in  the  formation  zone 
T  b  the  unbumed  temperature. 


Nagle 

oxidation: 


and  Strickland  [44]  proposed  the  following  relation  for  the  rate  of  soot 


ds 


(4.16) 


where: 

Pff^  b  the  partial  pressure  of  oxygen 

p,  b  the  density  of  the  soot  particle 

d,  b  the  diameter  of  the  soot  particle 

f(T,  Pf,  )  b  a  complex  Arrhenius  type  expression. 


Lee  [45]  proposed  the  following  soot  oxidation  model: 

—  -  d/  r  exp( -20,000 /J) 


(4.17) 
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Fraimorc  and  Jonc’s  [4d]  snsfcatcd  Aat  tkc  rate  oHxaoot  oiidatioa  is  gtvcn  by: 

^  -  -Pl^  P^o  exp(-19000/r)  (4.H) 

dt 


Hiroyasu  and  Kadota  [47]  proposed  the  following  equations  for  soot  formation  and 
oxidation: 

-  kl^-E^RT) 

p 

~  P'*  exp( -EJRT) 

_  dm^  dm^ 

~  dt  dt 

where: 

is  the  mass  of  formed  soot 
is  the  mass  of  oxidixed  soot 
is  the  mass  of  fhel  vapor 
mass  of  soot. 

E^  and  E^  are  assumed  to  be  1.25£-H)4  Kcal/Kmol  and  1.4E-H)4  Kcal/Kmol 
respectively.  Af  and  .4^  are  constants  which  are  determined  so  as  to  match  the  computed 
soot  with  the  experimental  results. 


dm^ 

“«F 

dm, 

dt 

dm^ 

~W 


Morel  [41]  recommended  empirical  correlations  for  soot  formation  and  oxidation 
directly  linked  to  the  combustion  modeL  The  soot  model  incorporates  only  those  variables 
calculated  by  the  combustion  model,  Le.,  temperature,  equivalence  ratio,  partial  pressure 
of  certain  species,  etc.  According  to  Mord  Ae  amount  of  soot  formation  is  described  by 


(4.1*) 


and  the  subsequent  soot  oxidation  in  the  burned  zone  is  described  by 

^  -  -«,»/(p,<y  aiK-«/r,)  j*“  ( 

where: 

5  is  the  mass  of  the  soot 

F  is  the  mole  fraction  of  available  oxygen  in  the  actively  burning  zone 
is  the  rate  of  fuel  burned  in  the  diffusion  burning  model 
A,  is  a  constant  »  038 
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Aj  is  a  constant  =  5000 
7)  b  the  temperature  in  the  formation  zone 
p,  is  the  soot  density  »  900  Kg/m’  (400  •  1500  Kg/m’) 

d,  is  the  diameter  of  the  elementaiy  soot  particle  entering  the  burned  zone  = 

0.012  -  0.032  pm 

is  the  partial  pressure  of  Oj  in  the  bnmed  zone 

Bf  is  a  constant  0.015 

Bj  is  a  constant  =  5000. 

Tesner  [3233]  suggested  a  kinetic  scheme  resulting  in  soot  formation  based  on  a 
chain*type  process  involving  radical  nuclei  from  which  soot  particles  will  later  grow. 

Magnussen  and  Hjertager  [28]  used  Tesner’s  soot  formation  model  to  analyze  a 
turbulent  combustion  flow,  and  added  a  model  of  soot  oxidation  which  considers  the 
behavior  of  soot  in  turbulent  flames.  Although  there  are  a  number  of  adjustable 
parameters  in  the  Tesner  model,  Magnussen  and  Hjertager  only  had  to  modify  Tesner*s 
coefficients  slightly  to  model  their  measurements  of  soot  formation  in  turbulent 
acetylene/air  flames.  These  models  are  described  in  more  detail  below. 

4.6.2  Soot  Model  Adapted  to  KIVA 

It  was  decided  to  use  the  soot  formation  model  developed  by  Tesner,  et  aL  [3233] 
and  the  sood  oxidation  model  developed  by  Magnussen  and  Hjertager  [28]  due  to  the  fact 
that  they  take  into  consideration  the  behavior  of  soot  in  turbulent  flames  and,  therefore, 
are  suitable  for  use  in  conjunction  with  the  eddy-break-up  combustion  model 

It  has  been  assumed,  following  Tesner,  that  soot  is  formed  from  a  gaseous  fuel  in 
two  stages,  where  the  first  stage  represents  formation  of  radical  nuclei,  and  the  second 
stage  represents  soot  particle  formation  from  these  nuclei  The  rate  of  formation  of 
radical  nuclei  is  expressed  by: 

+  (f-g)  n  -  .  H  .  N  (part/m’/s)  (^-21) 

where: 

/  is  a  linear  branching  coefficient 
g  is  a  linear  termination  coefficient 
g„  is  the  coefficient  of  linear  termination  of  soot  particles 
n  is  the  concentration  of  radical  nuclei  (part/m’) 

N  is  the  concentration  of  soot  particles  (part/m’) 

n„  is  the  rate  of  spontaneous  formation  of  radical  nuclei 

n„  b  expressed  by: 

(part/m’/s)  (^*22) 
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where: 


a,  is  a  constuit 

Cf  is  the  mass  concentratioii  of  fiiel  (Hx/m^ 

E  is  the  activatioB  temperature  (K) 

T  is  the  absolute  temperature  (K) 

R  is  the  universal  gas  constant. 

The  rate  of  soot  particle  formatton  is  expressed  by: 


(  a  -  b  N  )  H  (kg/m’/s) 


where: 

a  4  h  are  constants 

m,  is  the  mass  of  the  soot  particle  (Kg/part).  It  is  calculated  assuming  soot 
particles  to  be  spherical  and  having  an  average  diameter  d,  200  A  and 
density  p,  »  2000  Kg/m’. 

In  this  model,  the  soot  formation  is  limited  to  regions  of  unbumt  lueL 


Soot  oxidation  occurs  when  soot  particles  come  into  contact  with  oxygen.  The  soot 
combustion  model  developed  by  Magnussen  and  Hjertager  takes  into  account  tihe 
behavior  of  soot  in  turbulent  flames. 

For  regions  where  the  local  mean  soot  concentration  is  low  compared  to  the  oxygen 
concentration,  the  rate  of  soot  combustion  can  be  expressed  as: 


R^^- A,  C. 


(kg/m’/s)  {424) 


where: 


A  is  a  constant 

C,  is  the  local  mean  soot  concentration  (kg/m’) 
k  is  the  turbulent  kinetic  energy  em’/s’ 

8  is  the  rate  of  dissipation  of  turbulent  kinetic  energy  em’/s’. 

In  regions  where  the  oxygen  concentration  is  low,  the  oxygen  will  limit  the  rate  of 
soot  combnstioB.  The  soot  must  also  compete  for  oxygen  with  the  nnbnmed  fueL  In  this 
case,  the  rata  of  soot  combustion  is  given  by: 


C,  r,  +  C,  r. 


(kg/m’/s)  (425) 


where: 


is  the  stoichiometric  oxygen  requirement  to  bum  1  kg  soot 
is  the  stoichiometric  oxygen  requirement  to  bum  1  kg  fuel 
is  the  local  time  mean  oxygen  concentration  (kg/m’) 

is  the  local  time  mean  fuel  concentration  (kg/m’). 
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The  equation  that  gives  the  lower  reaction  rate  determines  the  local  rate  of  soot 
combustion.  Thus,  combining  equations  (4.24)  and  (4.25)  leads  to  the  expression: 


(kg/mVs)  i4M) 


The  local  number  of  radical  nuclei  can  be  assumed  to  be  reduced  by  combustion 
according  to: 


(part/m  Vs) 


(4.27) 


The  concentration  rate  of  soot  particles  is  given  by  the  difference  of  the  soot 
formation  rate  and  the  oxidation  rate,  i.e. 


^  -R 


(part/m  Vs) 


(4J«) 


Similarly,  The  concentration  rate  of  radical  nuclei  is  given  by: 


dll  Q  Q 

dt  ^  * 


(part/m  Vs) 


{4J9) 


The  simultaneous  solution  of  the  above  two  differential  equations  is  obtained  using 
a  fourth  order  Runge-Kutta  method. 

4.6.3  Soot  Model  Parameters 

Referring  to  studies  on  soot  formation  in  atmospheric  flames  and  shock  tubes 
showed  that  the  conditions  under  which  soot  formation  occurs,  as  well  as  the  quantity  of 
soot  formed,  vaiy  little  for  most  hydrocarbons.  Therefore,  the  constants  appearing  in  the 
soot  formation  equation  recommended  by  Tesner  et  al.  for  acetylene-oxygen  flames  have 
been  used  in  this  investigation  along  with  the  constants  recommended  by  Magnussen  [28] 
for  the  soot  oxidation  equations.  The  values  of  these  constants  are  given  below: 

a  =  1.6e5  /s 

f-g  =  100  /s 

g,  =  l.Oe-15  mVpart/s 

b  =  8e-14  mVpart/s 

Sg  =  7.254e30  part/kg/s 
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fyR>9«4  K 

p,  ■  2  g/cm^ 

d,  (the  soot  partkte  mean  diameter)  »  200  A 
The  local  mean  soot  concentration,  C.  (kg/m^,  is  defined  as: 

C,  *  p,  (kg/m’)  vol,  (m’/part)  N  (part/m’) 

where: 

p,  is  the  density  of  the  soot  particle  (kg/m’) 

vol,  is  the  volume  of  a  soot  particle  (m’/part) 

N  is  the  concentration  of  soot  particles  (part/m’). 

The  mass  of  a  soot  particle,  mp  (kg/part),  is  defined  as: 

mp  =  p,  (kg/m’)  vol,  (m’/part) 

»  2000  (kg/m’)  (4/3)n  (djlf  (m’/part) 

»:8J7758e-21  (kg/part) 

Therefore,  the  local  mean  soot  concentration,  C,  (kg/m’),  is  given  by: 

C,  =  8J7758e-21  (kg/part)  N  (part/m’) 

The  stoichiometric  oxygen  requirement  to  bum  1  kg  of  soot,  is  calculated  from 
the  foUowing  chemical  reaction  [48]: 

Soot  +  CO^ 


Therefore,  r,  is  given  by: 

r,  -  32/12  -  2.66667 

The  soot  modd  caicnlates  the  soot  concentration  in  (kg/m’).  The  calculation  stops 
at  90  degree  crank  ■«gl*  ATDC  where  the  calculated  average  cylinder  gas  temperature 
is  990  K,  and  the  average  cylinder  pressure  is  4.98  atm.  For  comparison  with  the 
experimental  data,  this  value  is  converted  to  exhaust  gas  conditions  using  the  following 
relation,  aMuming  ideal  gas  conditions. 

C,  (at  STP)  =  C,  (at  90**  ATDC)  *  ~p^  (kg/m’)  (4J0) 

Appendix  F  lists  the  soot  model  FORTRAN  source  code,  and  Appendix  G  shows  the 
modifications  made  to  the  KIVA-II  subroutine  NEWCYC,  and  subroutine  GLOBAL 
associated  with  the  soot  modeL 
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4.6.4  Results  of  Soot  Modeling 


Figures  4.17  and  4.18  show  the  computed  soot  formation,  soot  oxidation,  and  the 
soot  concentration,  which  is  the  different  between  the  formation  and  oxidation,  for 
hexadecane  and  dodecane  at  1500  RPM,  18  DBTDC  and  a  load  of  3.91  bar  BMEP. 
Figures  4.19  and  4.20  show  the  computed  soot  concentration  versus  crank  angle  for 
hexadecane  and  dodecane  at  three  different  BMEP  of  2.24,  3J5,  and  3.91  bar  at  1500 
RPM  and  an  injection  timing  of  18  DBTDC.  The  model  produces  the  correct  trend  of 
increasing  soot  concentration  with  increasing  load. 

The  rate  of  soot  concentration  versus  crank  angle  can  be  divided  into  three  regions. 
First,  soot  begins  to  form  in  the  fuel  rich  mixture  at  the  beginning  of  combustion  and 
during  the  injection  process  where  the  fuel-air  ratio  at  the  inner  zone  of  the  spray  is  high. 
Then  the  compression  stroke  continues  until  the  piston  reaches  TDC,  and  the  expansion 
stroke  begins.  During  this  period,  the  air  motion  (swirl)  is  low  and  the  rate  of  soot 
formation  is  higher  than  the  rate  of  soot  oxidation.  This  continues  to  about  20  degrees 
ATDC  and  the  soot  concentration  rises  sharply.  At  20  degrees  ATDC,  the  air  motion 
continues  increasing  at  a  higher  rate  as  the  piston  continues  moving  to  about  40  degrees 
ATDC.  This  results  in  increasing  the  rate  of  soot  oxidation.  It  is  observed  that  the 
maximum  soot  concentration  occurs  at  about  40  degrees  crank  angle  ATDC.  After  40 
degrees  crank  angle  ATDC,  the  air  motion  continues  to  increase  at  a  higher  rate  and  the 
small  amount  of  remaining  fuel  continues  to  bum.  On  the  other  hand,  the  soot  formed 
earlier  comes  in  contact  with  oxygen,  resulting  in  a  high  oxidation  rate,  and,  as  seen  from 
the  figures,  the  rate  of  soot  concentration  continues  to  decrease  to  about  60  degrees 
ATDC.  Finally,  the  mixture  reaches  equilibrium  and  the  amount  of  soot  concentration 
remains  constant  until  the  end  of  the  expansion  stroke. 

Figures  4.21  and  4.22  compare  the  experimental  soot  emission  for  the  baseline 
engine  with  the  computed  results.  The  Figures  show  the  effect  of  load  on  the  soot  emission 
at  STP  condition  (J=^9  K,  P^^l  atm)  for  hexadecane  and  dodecane  respectively.  The  soot 
model  produces  the  correct  trend  of  increasing  soot  concentration  with  increasing  load. 
As  shown  in  the  figures,  the  model  gives  fair  agreement  with  the  measured  soot  values. 
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SOOT  EMISSION  (G/M**3)  SOOT  EMISSION  (G/M*»3) 


Figure  4.17 

Soot  Fornation  aad  OiMatioa  Versui  Craak  Aagla  for  Haiadocaao 
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Figure  4.18 

Soot  Formatioii  and  Oiidation  Vcnui  Crank  Angia  for  Dodecano 
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Figure  4.19 

Soot  ConccntratioM  Vertuc  Creak  Aagie  hr  Hexadecaae 
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Figure  4.20 

Soot  Conceatration  Venua  Creak  Aaglo  for  Dodccaae 
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EMISSION  (G/M**3)  SOOT 


BMEP  (BAR) 

Figure  4Jlt 

Soot  Emiuioa  Versus  Load  for  Houdccano 


BMEP  (BAR) 

Figure  4X1 

Soot  Emissiou  Versus  Load  for  Dodccauo 
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4.7  THERMAL  MODELING  OF  CERAMIC  COATINGS 


Two  of  the  input  parameters  required  for  KIVA-D  are  die  surface  temperature  of 
the  combustion  chamber  walls  and  the  intake  air  temperature.  For  the  baseliae  case,  the 
temperature  of  all  surfaras  was  assumed  constant  and  uniform  at  425  K,  and  the  intake 
air  temperature  was  assumed  to  be  315  K.  The  effect  of  applying  a  thin  ceramic  coating 
to  combustion  chamber  surface  is  to  increase  both  the  average  surface  temperature  and 
the  variation  in  surface  temperature  relative  to  the  baseline  metallic  surface.  To  model 
this  effect,  the  coated  surfaces  were  assumed  to  have  a  cyclic  temperature  profile,  as 
depicted  in  Figure  4.23,  consisting  of  a  100  K  variation  centered  on  an  average  value  of 
475  K.  This  temperature  profile  is  based  on  results  of  experimental  and  analytical  studies 
related  to  ceramic  coated  surface  temperatures  as  cited  in  references  3,  6,  and  7.  The 
uncoated  surfaces  of  the  combustion  chamber  were  kept  at  the  same  value  as  for  the 
baseline  model,  i.e.  425  K,  following  a  study  by  Smavik  {8)  who  observed  ne^igible 
temperature  rise  in  the  uninsulated  components.  The  intake  air  temperature  at  the  start 
of  compression  was  also  increased  to  account  for  both  the  mixing  with  higher  exhaust  gas 
temperature  residuals  caused  by  insulated  surfaces  and  heat  transfer  from  the  higher 
temperature  surfaces  to  the  intake  air.  For  the  coated-piston  case  and  the  all>surface« 
coated  case,  325  K  was  assumed  for  the  intake  air,  while  for  the  coated-head  case,  320  K 
was  assumed.  These  values  were  based  on  the  experimental  volumetric  efficiency  data. 


4.8  THERMAL  MODELING  RESULTS 


KIVA  results  are  presented  in  Figures  4.24-4.27  for  four  different  sets  of  operating 
conditions.  All  results  shown  are  for  hexadecane  fuel  at  1500  RPM.  Figure  4.24  is  for 
a  load  of  8  NM  (223  KPa  BMEP)  and  a  timing  of  22  DBTDC.  Figure  4.25  is  for  a  load 
of  12  NM  (335  KPa  BMEP)  and  a  timing  of  20  DBTDC.  Figure  4.26  is  for  a  load  of  12 
NM  (335  KPa  BMEP)  and  a  timing  of  22  DBTDC.  Figure  4.27  is  for  a  load  of  14  NM 
(391  KPa  BMEP)  and  a  timing  of  22  DBTDC.  Each  figure  includes  plots  of  cylinder 
pressure,  cylinder  gas  temperature,  mass  of  burnt  fuel  and  NO  concentration  for  the 
specified  load  and  injection  timing.  Each  plot  compares  the  three  ceramic  coated  cases 
with  the  baseline  case.  The  effect  of  the  various  insulation  schemes  on  pressure, 
temperature,  rate  of  fuel  combustion  and  NO  emission  is  discussed  in  the  following 
sections. 
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4.8.1  Cy Under  Gas  Pressure 


The  P*6  diagrams  in  Figures  4.24(a)-  4.27(a)  are  simulations  from  KIVA  for  the 
baseline  case  and  the  three  insulated  engines.  In  each  case  the  peak  cylinder  pressure  is 
highest  for  the  baseline  engine,  followed  by  the  coated-head  case,  the  all-surface-coated 
case  and  the  coated-piston  case.  The  coated-piston  case  always  produces  the  lowest  peak 
cylinder  pressure.  These  KIVA  results  agree  weU  with  the  experimental  P-6  diagrams 
shown  in  Figure  3.1.  In  general,  the  peak  cylinder  pressure  is  a  function  of  the  ratio  t^  ft, 
(premixed  combustion  duration/diffusion  combustion  duration).  The  higher  the  ratio,  the 
higher  the  peak  cylinder  pressure,  and  vice  versa.  This  leads  to  the  conclusion  that 
selective  insulation  of  the  combustion  chamber  walls  can  modify  the  combustion 
characteristics  of  the  engine,  causing  the  peak  cylinder  pressure  to  vary  according^.  It 
also  impUes  that  the  ignition  delay  for  the  coated-piston  engine  is  the  shortest  and  the 
ignition  delay  for  the  baseline  engine  is  the  longest  of  the  four  engine  builds  investigated. 
This  is  substantiated  by  the  experimental  data  in  some  cases  but  not  in  others,  as 
discussed  in  Section  3.5. 

4.8.2  Cylinder  Gas  Temperature 

KIVA-n  calculates  the  gas  temperature  in  each  computational  cell  at  each  time 
step.  Relations  were  added  in  this  study  to  further  calculate  the  mass  average 
temperature  of  the  entire  cylinder  contents  based  on  an  energy  balance  at  each  time  step. 
Figures  4.24(b)  -  4.27(b)  compare  the  calculated  average  gas  temperature  versus  crank 
angle  for  the  baseline  and  three  selective  insulation  schemes.  The  following  general 
observations  can  be  made  from  these  figures:  1)  During  compression,  the  cylinder  gas 
temperature  u  higher  for  the  coated  cases  than  for  the  baseline,  2)  the  gas  temperature 
for  the  various  ceramic  coated  cases  begins  to  rise  earlier  but  at  a  slower  rate  than  the 
baseline  case,  3)  the  maximum  average  cylinder  gas  temperature  is  lowest  for  the  coated- 
piston  case  and  highest  for  the  baseline  case,  with  the  coated-head  and  totally  insulated 
cases  between. 

These  observations  may  be  attributed  to  tbe  effect  of  the  selective  coating  schemes 
on  the  ignition  delay  period.  The  initial  higher  gas  temperature  is  due  to  the  higher 
surface  temperature  of  the  ceramic  surfaces.  The  eariier  rising  temperature  is  indicative 
of  a  shorter  ignition  delay  which  results  in  a  smaller  premixed  combustion  fraction  and 
a  longer  diffusion  combustion  period  resulting  in  a  lower  average  peak  cylinder 
temperature  and  a  higher  gas  temperature  during  the  expansion  stroke.  The  results 
discussed  above  for  peak  cylinder  pressure  are  also  in  support  of  this  reasoning. 
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4.8J  Mass  of  Burnt  Fuel 


Figures  4.24(c)  -  4.27(c)  depict  the  mass  of  fuel  hir  mt  during  the  combustion 
process.  In  general,  it  can  be  noted  that  the  rapid  rate  of  co.iibustion  associated  with  the 
premixed  combustion  phase  occurs  earlier  and  consumes  less  fuel  for  all  of  the  coated 
cases  relative  to  the  baseline  case,  indicative  of  shortened  ignition  delay  and  decreased 
t^t^.  This  is  in  keeping  with  the  lower  cylinder  pressure  and  temperature  predictions. 
Typically  it  is  the  coated-piston  case  that  displays  the  lowest  rate  of  combustion. 

4.8.4  Nitric  Oxide  Emission 

The  rate  of  NO  formation  during  combustion  for  the  baseline  and  three  selective 
insulation  schemes  is  shown  in  Figures  4.24(d)  -  4.27(d).  It  is  noted  that  the  coated- 
piston  case  typically  produces  the  lowest  NO  emission  consistent  with  the  lowest  rate  of 
combustion  and  cylinder  temperature.  The  baseline  engine  typically  produces  the  highest 
NO  emission  except  for  the  case  shown  in  Figure  4.27(d)  for  14  NM  (391  KPa  BMEP) 
load  and  22  DBTDC.  In  this  case  the  coated  head  produced  a  higher  NO  value  than  did 
the  baseline.  The  experimental  values  of  NO  for  these  cases  are  indicated  by  the  data 
points  shown  along  the  right  side  Y-axis.  Reasonable  agreement  between  experimental 
and  predicted  values  is  shown  in  most  cases  and  the  correct  relative  values  of  the  coated 
cases  versus  baseline  are  obtained  except  for  the  low  load  condition  shown  in  Figure 
4.24(d)  where  agreement  is  not  obtained  for  the  insulated  bead  case.  Figures  4.28  and 
4.29  further  illustrate  reasonable  agreement  between  experimental  values  of  NO  and  those 
predicted  with  KIVA,  especially  at  higher  values  of  load.  At  low  loads  KIVA 
underpredicted  NO  by  about  25  percent. 
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KIVA  Results  for  1500  RPM,  22  DBTDC  and  223  KPa  BMEP  Load  (8  NM  Torque) 
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Figure  4J5(a) 
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Figure  4.25 

KIVA  Results  for  1500  RPM,  20  DBTDC  and  335  KPa  BMEP  Load  (12  NM  Torque) 
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Figure  4.26 

KIVA  Results  for  1500  RPM,  22  DBTDC  and  335  KPa  BMEP  Load  (12  NM  Torque) 
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Moaa  of  Burnt  Fuol  (G) 


Figure  4J7(a) 

Cylinder  Gas  Pressure  Versus  Crank  Angle 
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Mass  of  Fuel  Burnt  Versus  Crank  Angle 


Crank  Angle  (OEG) 


Figure  4J7(b) 

Cylinder  Gas  Tanigaratam  Versos  Crank  Angle 


Figure  417(d) 

Nitric  Glide  Concentration  Verms  Crank  Angle 


Figure  4.27 


KTVA  Results  for  1500  RPM,  22  DBTDC  and  391  KPa  BMEP  Load  (14  NM  Torque) 
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Figure  D1 

Soot  Conccntnition  Versus  Load  for  Hexadcoine  at  16  DBTDC 
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Figure  D3 

Soot  Concentration  Versus  Load  for  Hezadecane  at  20  DBTDC 
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Figure  D4 

Soot  Concentration  Versus  Load  for  Hexadecane  at  22  DBTDC 
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Figure  D5 


Soot  Concentration  Versus  Load  for  Dodcauie  at  16  DBTDC 


IroHon  (8/m*3)  Sool  ConcantroHon  (g/m*3) 


1000  RPM  &  18  OeTDC 

OOOECAIC 


1500  RPM  ft  18  OBTDC 


300  400  SOO 

200  300  400 

OMP  (XPa) 

BtKPlKPo) 

Fig  D6(ft) 

FigD6(b) 

2000  RPM  ft  18  OBTDC 

2500  RPM  ft  18  OBTDC 

I 


1000  RPM  ft  20  OeTDC 


1500  RPM  ft  20  OeiOC 


500  400 

eiO>(KI*a) 


^<■2  Skim. 
•o  -«• 

V  OMMIM 

^  ctMr.  «  i 


500  400 

9KP  (KPo) 


Fig  D7(ft) 


Fig  D7(b) 


2000  RPM  ft  20  OBTOC 


2500  RPM  ft  20  DBTOC 


^  . 


500  400  SOO 

SMF(KI>a) 


Fig  D7(c) 


Fig  D7(d) 


Figure  D7 


Soot  Concentration  Versus  Load  for  Dodecane  at  20  DBTDC 
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Figure  D8 


Soot  Concentration  Versus  Load  for  Dodecane  at  22  DBTDC 
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Figure  El 

Exhaust  Gas  Temperature  Versus  Load  for  Hex;adecane  at  16  DBTDC 
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Figure  E2 

Exhaust  Gas  Temperature  Versus  Load  for  H^decane  at  18  DBTDC 
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Figure  E3 

Exhaust  Gas  Temperature  Versus  Load  for  Hmuidecane  at  20  DBTDC 
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Figure  E5 

Exhaust  Gas  Temperature  Versus  Load  for  Podecane  at  16  DBTDC 
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Figure  E8 

Exhaust  Gas  Temperature  Versus  Load  for  Dpdccane  at  22  DBTDC 
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